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Abstract 

We discuss the interpretation of present and possible future experimental constraints on 
B s — > fi + fi~ decay in the context of the constrained minimal extension of the Standard 
Model (CMSSM) with universal scalar masses. We emphasize the importance of including 
theoretical and other experimental uncertainties in calculating the likelihood function, which 
can affect significantly the inferred 95 % confidence-level limit on the CMSSM parameters. 
The principal uncertainties are the B s meson decay constant, m t and m&. The latter induce 
uncertainties in the mass of the neutral Higgs boson that dominates the B s — > fi + decay 
amplitude at large tan /3, reducing the CMSSM region excluded by present and possible 
future limits from the Fermilab Tevatron collider and the LHC. 
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1 Introduction 



The prospects for new physics searches at the LHC and other future colliders are already 
constrained by rare processes that are sensitive to small deviations from the Standard Model. 
A prime example of such a low-energy constraint is b — > 57 decay [1,2]. This, the anomalous 
magnetic moment of the muon, — 1 [3], and the Higgs mass, [4], are among the most 
important indirect constraints on extensions of the Standard Model, such as the minimal 
super symmetric extension of the Standard Model (MSSM). 

The decay B s — > fi + fi~ has been emerging as another interesting potential constraint 
on the parameter space of models for physics beyond the Standard Model, such as the 
MSSM [5-8]. The Fermilab Tevatron collider already has an interesting upper limit on 
the branching ratio for B s — > [9], and is expected soon to increase significantly its 

sensitivity to B s — > fi + fi~ decay, whilst the LHC sensitivity will reach down to the Standard 
Model rate [10]. Since its branching ratio is known to increase rapidly for large values of the 
ratio of Higgs v.e.v.'s, tan/3 [11-13], increasing like its sixth power, these present and future 
sensitivities are particularly important for models with large tan (3. 

In view of its potential importance for the MSSM, it is important to treat the B s — > fi + {i~ 
decay constraint with some care, as has already been done for the b — > 57, — 2 and rrih 
constraints. In each of these cases, the interpretation depends on uncertainties in theoretical 
calculations, which should be propagated carefully and combined with the experimental 
errors if limits are to be calculated at some well-defined confidence level, or if a global fit to 
MSSM parameters is to be attempted. 

Both of these issues are important also in the case of B s — > fi + fi~ decay. As concerns 
the theoretical uncertainties, it is important to include consistently all the available one-loop 
MSSM contributions, and avoid any approximate treatments of any individual terms, since 
cancellations are potentially important, and also to include knowledge of the higher-order 
QCD corrections to the lowest-order MSSM loop diagrams. As we discuss in this paper, 
other sources of error and uncertainty are also important. These include, for example, the 
uncertainties in the B s meson parameters, principally the decay constant fs s - Since the 
B s — > fi + fi~ decay rate is dominated by the exchange of the pseudoscalar Higgs boson A, 
the value of itia is also an important uncertainty. 

Bounds on supersymmetry are often explored in a constrained model with universal soft 
supersymmetry-breaking scalar masses mo, gaugino masses m^j and trilinear couplings Aq 
at some GUT input scale. In this CMSSM, the values of tua and magnitude of the Higgs 
mixing parameter jj, (but not its sign) are in principle fixed by the electroweak vacuum 



1 



conditions. However, these predictions are necessarily approximate. For example, the value 
of rtiA predicted as a function of the independent parameters mi/ 2 ,m , A and tan/3 has 
significant uncertainties associated with the lack of precision with which the heavy quark 
masses m t and m b are known, as we discuss extensively later in this paper. Moreover, rather 
different values of m A would be predicted in models where the universality assumptions of 
the CMSSM are relaxed. For example, much smaller values of m A are attainable in models 
with non- universal Higgs masses (NUHM). 

When interpreting experimental upper bounds (or measurements) within any specific 
model such as the CMSSM, care must be taken to incorporate the uncertainties in auxiliary 
parameters such as fs s , m t and mj. These should be propagated and combined with the 
experimental likelihood function when quoting sensitivities in, e.g., the (mi/ 2 ,mo) plane at 
some fixed level of confidence. Moreover, one must also be aware of model dependences 
within the assumed framework such as the value of A in the CMSSM, as well as the effects 
of possible deviations from the model framework such as non-universal Higgs masses. 

We exemplify these points in a discussion of uncertainties in the interpretation of the 
present and likely future sensitivities of the Fermilab Tevatron collider and the LHC to 
B s — > decay, assuming /i > as preferred by — 2. In particular, we show that the 

uncertainties in f Bs , m t and m b each shrink significantly the regions that might otherwise 
appear to be excluded by the present limit, or might appear to be if the decay is not discovered 
at the likely future sensitivity. We compare the resulting B s — > constraints with other 

existing constraints such as b — > sj, discussing how they vary with A and commenting on 
the situation within the NUHM. 



2 Calculation of B s — ► /x + fi Decay 

The branching ratio for the decay B s — > ji + /i~ is given by 
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where the one-loop corrected Wilson coefficients C s ,p are taken from [13] and Ca is defined 
in terms of Y(x t ), following [14], as C A = Y(x t )/ sin 2 9 W where 
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The function Y(x t ) incorporates both leading [15] and next-to-leading order [10] QCD cor- 
rections, and m t {m t ) is the running top-quark mass in the MS scheme. The precise value 
of m t (m t ) depends somewhat on the set of supersymmetric parameters and our choice of 
the physical top quark mass m t = 178 ± 4 GeV [16] that we use in this paper. The Wilson 
coefficients Cs,p have been multiplied by 1/(1 + e b ) 2 , where e b incorporates the full super- 
symmetric one-loop correction to the bottom-quark Yukawa coupling [17-19]. It is known 
that, since e b is proportional to tan/3, this correction may be significant in the large-tan (3 
regime we study here [5,6]. These corrections to e b are flavour independent. In addition, it is 
important to include the flavour-violating contributions arising from the Higgs and chargino 
couplings at the one-loop level. These corrections result in effective one-loop corrected values 
for the Kobayashi-Maskawa (KM) matrix elements [12,20], which we include as described 
in [21,22]. In particular, these corrections modify the Wilson coefficients involved in Eq. (1), 
as can be seen in Eqs. (6.35) and (6.36) in [21] or in Eq. (14) in [22]. In addition, the 
latter study includes the effects of squark mixing, which are included here as well. Below, 
we discuss the behaviour of the dominant contribution to the Wilson coefficients, mainly to 
illustrate their dependence on the pseudo-scalar Higgs boson mass. 

The Wilson coefficients Cs,p receive four contributions in the context of MSSM, due 
to Higgs bosons, counter-terms, box and penguin diagrams. The Higgs-boson corrections 
were calculated in [14], and the rest of the supersymmetric corrections in [11,12]. The 
full one-loop corrections have been studied and presented comprehensively in [13]. In our 
numerical analysis, we implement the full one-loop corrections taken from this work as well 
as the dominant NLO QCD corrections studied in [23], as well as the flavour- changing gluino 
contribution [22,24] . The Higgs-boson, box and penguin corrections to Cs,p are proportional 
to tan 2 (3, while the counter-term corrections dominate in the large-tan f3 limit, as they are 
proportional to tan 3 (3. 

In order to understand the behaviour of the branching ratio in the (ffii/2,mo) plane in 
the context of the CMSSM, we focus attention on the counter-terms which are mediated by 
A,H,h exchange as seen in Eqs. (5.1) and (5.2) of [13]. As seen in Eq. (5.13) of [13], the 
B s — > fi + fi~ decay amplitude oc l/m\ in the large-tan (3 limit, and the decay rate oc l/m\. 
This underlines the importance of knowing or calculating as accurately as possible. 

The counter-term contribution to Cs,p is given by [13] 

y/zm w m A . =1 a=1 =1 



3 



where 

V t mn = ~ ,= 1 2a [V2M w V tl (T u -% a - (MuUV^F^UX^, (4) 
2v 2 sin O w 

and My = diag(m n , m c , m t ). U and V are the chargino diagonalization matrices, Y Ul and 
Y Ur are 6x3 squark diagonalization matrices, and D 3 (x) = x\nx/(l— x). Additionally, y a i is 
defined in Eq. (5.10) of [13] as y ai = m 2 -Jm 2 x +, where m\ a = {m\ L ,m\ h ,m\,m\ R ,m\ R ,m\^ . 
Finally, \ mn = V mb V* s /V tb V*. 

We can split the counter-term contribution into two terms: one that is proportional 
to My/ and another that is proportional to m t . Starting with the term whose numerator 
depends on My/, it is easy to see that the non- vanishing terms stem from the following 
combinations of indices: {a,n,m} = {111,222,333,633} and i = 1,2. However, the first 
term {111} is suppressed, as it is proportional to An = V ub V us /Vtb Vt s — —0.022, whereas 
the second is not suppressed, because it is proportional to A22 = V cb V cs /Vtb Vt s — — 1. Nor 
are the third and fourth terms suppressed, as they are multiplied by A33 = 1. Thus, the part 
of the counter-term contribution to the Wilson coefficient that is oc M w is 

c CT,Mw = f y u jj^ ^ D ^ y2i) + ( CQS 2 iD ^ yzl ) + sin 2 Q. D 3 (y 61 ))] 

V21 U 22 [A 22 £3(2/22) + A 33 (cos 2 ^£3(2/32) + sin 2 ^£3(2/62))] } , (5) 



+ m x 2 



where 



_ rrif, tan 3 p 
/ "4M^sin 2 ^m 2 4 ' [) 
and we have ignored in (5) terms that are proportional to An. The unitarity of the KM 
matrix implies that An + A 2 2 + A 33 = 0, which for small An means A 2 2 = — A 33 , resulting in 
the suppression of C° s ™ w . 

Turning now to the terms that increase with the charge-2/3 quark masses, we see that the 
terms with n = 3 (the top-quark contributions) dominate the first- and second-generation 
terms in T°; mn . Specifically, the dominant terms have {a,n, m} = {333,633}. In addition, 
we notice that the i — 2 part is important, since it is proportional to V22 ^22, while the i — 1 
term is proportional to V\ 2 U± 2 . Hence it is sufficient to take 

C C s T ' mt = m t f (^) m xt [D(y 32 ) - £(2/62)] , (7) 

where we set V 2 2 — U 22 — 1- The expression (7) is the approximation derived in [6]. The 
GIM cancellation in (5) means that the counterterm contribution to the Wilson coefficients, 
which is the dominant one, can be approximated relatively well by (7). However, in our 
analysis we use the full expressions given in [13], as well as the gluino correction and the 
flavour-violating corrections to Kobayashi-Maskawa matrix elements as described above. 
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3 Dependence of tua on m t and on m& 



It is clear from the discussion above that the mass of the pseudoscalar Higgs boson, m^, is 
an important ingredient in calculating the branching ratio for the decay B s — > fi + fi~ , since 
it enters in the fourth power in the denominators of the Wilson coefficients C,s and Cp in 
(1). Therefore, to further our discussion of the uncertainties in the B s — > fi + fi~ branching 
ratio, we now discuss the uncertainties in the calculated value of uia- 

The electroweak symmetry breaking conditions may be written in the form: 

m \ = m k + m k + V + A A (8) 

and 

2 _ rn? Hl - m 2 H2 tan 2 (3 + \m\{\ - tan 2 g) + Ag 
" " tan 2 /5-l + A( 2 ) ' W 

where A^ and A^ 1,2 -* are loop corrections [19,25-28]. The exact forms of the radiative correc- 
tions to fj, and rriA are not needed here, but it is important to note that, at large tan/3, the 
dominant contribution to contains a term which is proportional to hf tan (3 2 , whereas the 
dominant contribution to m 2 A contains terms proportional to h 2 tan (3 and h 2 tan (3. There- 
fore, the m 2 H2 term along with a piece proportional to h 2 in A^ are dominant in (9) in the 
large-tan (3 regime, so \i depends rather mildly on m&. 

We illustrate in Fig. 1 the logarithmic sensitivity of tha, namely Am^/m^, to m& and 
m t along slices through the (mi/ 2 ,m ) plane for tan/5 = 57, A = and /i > 0. We use 
as representative errors Am t = 1 GeV and Am b = 0.1 GeV. Panel (a) shows the effect of 
varying m for fixed mi/ 2 = 300 GeV, and panel (b) shows the effect of varying mi/ 2 for 
fixed m = 400 GeV. 

One can understand the behaviour depicted by employing (8), (9) and the renormalization- 
group equations (RGEs) that govern the evolution of the m Hi from M GUT to M z . The 
one-loop RGE for m 2 Hl depends on h\>: 

° m * - 1 6g 2 2 M| - %g\ M 2 ) + 6|/^| 2 (m 2 , + m 2 + m 2 + A 2 )} , (10) 



shig " i6^ I™ J " 2 ~ 3 yi ^ U|,tb| ^ " tD3 ^ 

where we ignore, for simplicity, the h T contribution, while the RGE for m 2 H2 depends on h t : 

f^f = 1^ {(-^ M l - \al Ml) + 6N VJ 3 + < + < + A 2 )} • (ii) 

We see from (10) that, as m b increases, the value of m 2 Hi at M z decreases, tending to 
decrease uia (8). More importantly, the term proportional to /i 2 in A^ in (8) is proportional 
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m (GeV) m 1/2 (GeV) 



Figure 1: The sensitivities of rriA to m t and m b , assuming Am t = 1 GeV and Arrib = 
0.1 GeV, along slices through the (mi/ 2 ,mo) CMSSM plane for A = and tan/3 = 57. 
Panel (a) fixes mu2 = 300 GeV and varies m , while panel (b) fixes m = 400 GeV and 
varies m-i/2- 

to m~ j (m~ —mf ) log(m~ /m~ )— 1, which is negative across the (mi/2, mo) plane for tan (3 = 
57 and A = 0, with a magnitude that decreases with mu2- Thus, as m& increases, we obtain 
a further decrease in m^. This term in fact provides most of the numerical dependence of 
rriA on m&. Since it is enhanced by tan/3, the m& dependence becomes milder for smaller 
values of tan /?. 

In addition, the effect on is augmented if the role of the bottom Yukawa coupling in 
the RGE (10) is enhanced, which occurs at large ttiq. This increases the sensitivity of 
to m&, as seen in Fig. 1(a). In contrast, when m.1/2 is increased, the sensitivity to is 
diminished and, at the same time, the gaugino part of the RGE is enhanced. Both effects 
lead to a reduced change in at large TO1/2, as can be seen in Fig. 1(b). 

Turning now to the dependence of tua on m t , we see that the evolution of m^ 2 shown in 
(11) is similar to (10), apart from the substitution of h t and analogous mass substitutions. 
As m t increases, the value of m^ 2 at Mz is driven to larger negative values. However, the 
change in tjia is dominated by the change in /x, which grows with an increase in m t . The net 
effect is an increase in m^, as seen in Fig. 1, which increases with m and decreases with 
rai/2, as seen in panels (a) and (b), respectively. 

Since the B s — > ji + n~ decay rate depends on the fourth power of m^, the sensitivity of 
rriA to both mt, and mt displayed in Fig. 1 can lead to a large uncertainty in B s — > fi + /i~ 
for tan/3 = 57. We have also evaluated the sensitivities to m t and mb for tan/3 = 40. These 
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sensitivities do not vary significantly with TO1/2 nor with m . They are always smaller than 
for tan/3 = 57, and the difference is rather substantial for large tuq and small ra^. 
Numerically, in the following analysis we assume 

Am t = 4 GeV, Am b = 0.11 GeV, (12) 

with central values for the physical pole mass m t = 178 GeV and the running mass m^ s {mb) = 
4.25 GeV. The first of the uncertainties in (12) is taken directly from measurements at the 
Fermilab Tevatron collider, and may be reduced by a factor of 2 to 4 by future measure- 
ments there and at the LHC. The following analysis shows that such reductions would be 
most welcome also in the analysis of B s — > ji + /i~ decay. The uncertainty in m& is taken from 
a recent review [29] that combines determinations from bb systems, 6-flavoured hadrons and 
high-energy processes. Our 1-a range given in [29] is contained within the preferred range 
quoted by the Particle Data Group [30] , and is very similar to the ranges quoted recently by 
the UKQCD group [31] in the unquenched approximation and in the review given by Rakow 
at the Lattice 2004 conference [32]. 

It is easy to see how important these uncertainties could be. For example, when Am A /m A = 
0.05 for Am t = 1 GeV, which occurs when tan/5 = 57 for (m 1/2 ,m ) = (300, 900) or (100, 
400) GeV, the change in m A for |Ara t | = 4 GeV is ±0.2m A , corresponding to a change in 
the A contribution to the B s — > fi + decay rate by a factor 2.07 or 0.41, depending on the 
sign of Am t . 

We note in passing that both the CDF [33] and DO [34] collaborations have recently pub- 
lished new upper limits on Higgs production at the Fermilab Tevatron collider. In particular, 
the DO limit [34] is relevant to the MSSM at very large tan/3 and small m A . However, such 
small values of m A are already excluded at large tan (3 in the CMSSM by other constraints 
such as b — > S7 and the lower limit on m h , so the present DO limit does not further restrict 
the part of the CMSSM parameter space of interest here. 

4 The Effects of Uncertainties on B s — ► /x + /x~ Limits 
in the CMSSM 

In order to assess how important these auxiliary uncertainties may be in the interpretation 
of B s — > ji + fi~ experiments, we display in Fig. 2 their individual effects on the present 
B s — > /i + iT constraint in the (■m l / 2 ,m ) plane of the CMSSM for A = 0, // > and 
tan p = 57. This is close to the largest value of tan (3 for which we find suitable electroweak 
vacua in generic domains of the (mi/ 2 ,mo) plane, and so maximizes the potential impact of 
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the B s — > constraint, which increases asymptotically as the sixth power of tan/5. The 

dark (brick) shaded regions in the bottom-right corners of each panel are excluded because 
there the lightest supersymmetric particle (LSP) would be the charged f\. The pale (blue) 
shaded strips are those favoured by WMAP, if all the cold dark matter is composed of LSPs. 
The supersymmetric spectrum and relic density calculations have been descibed elsewhere 
(see e.g. [35]). The near-vertical dashed (black) lines at small mi/ 2 are the constraint imposed 
by the non-observation of charginos at LEP, and the near-vertical dash-dotted (red) lines 
are those imposed by the non-observation of the lightest MSSM Higgs boson, as calculated 
using the FeynHiggs code [36]. 

The medium (green) shaded regions are excluded by the rare decay b — > sj. The branch- 
ing ratio for this has been measured by the CLEO, BELLE and BaBar collaborations [1]. 
The theoretical prediction of b — > 37 [2,21,37] contains uncertainties which stem from the 
uncertainties in m&, a s , the measurement of the semileptonic branching ratio of the B meson, 
and the effect of the scale dependence. In particular, the scale dependence of the theoretical 
prediction arises from the dependence on three scales: the scale where the QCD corrections to 
the semileptonic decay are calculated and the high and low energy scales relevant to b — > 57 
decay [38]. These sources of uncertainty can be combined to determine a total theoretical 
uncertainty 1 . The experimental measurement is converted into a Gaussian likelihood and 
convolved with a theoretical likelihood to determine the total likelihood [40], which is used 
to calculate the excluded region at 95% CL. It is important to note that the dependence of 
this excluded region on m^mi,, and m t is quite weak in comparison, as we have checked 
numerically. 

Finally, the ellipsoidal contours represent the nominal B s — > ji + fj,~ branching ratio, 
calculated (like all the others) using the current central values m t = 178 GeV and mff s (m b ) = 
4.25 GeV. The numerical labels for the two outer solid lines are exponents in the branching 
ratio: 1CT 7 , 1CT 8 , the thinner-dashed line is for 2 x 10~ 8 , the thicker-dashed line for 5 x 1CT 8 . 

The most stringent experimental upper bound on the B s — > /i + /i~ branching ratio is 
that given by an updated CDF measurement: 1.5 x 1(T 7 (2.0 x 10~ 7 ) at the 90% (95%) 
CL [9] . The innermost thick solid line of Fig. 2 is the contour for the present nominal 95% 
CL experimental upper limit of 2.0 x 10~ 7 . 

Panel (a) of Fig. 2 displays as a blue shaded region the effect on the interpretation of the 

present experimental limit of including the uncertainties in the B s meson parameters. The 

1 According to a recent analysis [39], these theoretical uncertainties may be significantly larger, resulting 
to a weaker bound on the masses of supersymmetric particles. 
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tan (3 = 57 , (i > tan [3 = 57, n>0 




m 1/2 (GeV) m 1/2 (GeV) 



Figure 2: The effects of auxiliary uncertainties on the region of the (mi/2, mo) plane for 
Aq = 0, fi > and tan/3 = 57 currently excluded by the Fermilab Tevatron collider, (a) 
The effect of B s meson uncertainties alone, principally that in /b s . (b) These uncertainties 
combined with the uncertainty Am< =4 GeV. (c) The B s meson uncertainties combined with 
the uncertainty Amj = 0.11 GeV. (d) All uncertainties combined. The various contours and 
shadings in the (m^jfflo) plane are explained in the text. 
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most important uncertainty is that in the decay constant fs 3 , for which we assume [41] 



f Ba = 230 ± 30 MeV. 



(13) 



In addition, we use [30] 



m Bs = 5369.6 ± 2.4 MeV, 



r Bs = (1.461 ± 0.057) x 10~ 12 s. 



(14) 



In calculating the uncertainty we add quadratically the uncertainties that result from these 
errors as well as those in the KM elements. We see that incorporating these uncertainties 
does not change the overall shape of the excluded region, but does shrink it slightly. There 
may be some possibility to reduce the uncertainty in f Bs in the foreseeable future, but in 
this analysis we retain it fixed in the following panels and other figures. 

The blue shaded region in panel (b) of Fig. 2 incorporates the present uncertainty in 
m t) assumed to be Am t = 4 GeV, which is propagated through the CMSSM calculation of 
rriA as discussed in the previous Section. We see that this uncertainty is more important 
for larger m , truncating the upper part of the exclusion domain. This effect can readily be 
understood from panel (a) of Fig. 1, where we saw that Am t has a particularly important 
effect on rriA at large m . 

The blue shaded region in panel (c) of Fig. 2 shows the parallel effect of the uncertainty 
in m b) assumed to be Am b = 0.11 GeV, which is also propagated through the CMSSM cal- 
culation of rriA as discussed in the previous Section. This uncertainty is also more important 
for larger m , providing an independent mechanism for truncating the upper part of the 
exclusion domain. This can also readily be understood from panel (a) of Fig. 1, where we 
saw that Amj also has a particularly important effect on rriA at large mo- 

There is also some tendency in both panels (b) and (c) of Fig. 2 for the exclusion domain 
to separate from the axis ~ 100 GeV, particularly at large m . This can be understood 
from panel (b) of Fig. 1, where we see that the effects of both Am t and Amj on are 
enhanced when mi/2 ^ 200 GeV. The uncertainties in each of m t and m b become particularly 
important when m 1 / 2 is small and m large, as seen separately in panels (b) and (c) of Fig. 2. 

The similar tendencies in panels (b) and (c) of Fig. 2 are reinforced when we combine the 
uncertainties in m t and m b , as shown by the blue shaded region in panel (d). We find that 
B s — > fi + /i~ decay is currently unable to exclude any value of m Q above about 350 GeV, 
whereas the exclusion region would have extended up to m ~ 450 GeV if the auxiliary 
uncertainties had not been taken into account, and ~ 400 GeV if either of the m t or mj 
uncertainties had been ignored. On the other hand, the reduction in the excluded range of 
my 2 at lower m is less important, typically < 30 GeV. 
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m 1/2 (GeV) m 1/2 (GeV) 



Figure 3: The disallowed regions in the (my2,mo) plane for fi > and tan/3 = 57, for (a) 
A = 2nii/2 and (b) A = — 2m 1 /2- The various contours and shadings in the (m^jMo) 
plane are as explained in the text describing Fig. 2. 

We observe that the region currently excluded by B s — > fi + [i~ is always included within 
the region already excluded by b — > S7 and/or m^, even without including the auxiliary 
uncertainties. The same is even more true for smaller values of tan/3: in the case tan/3 = 40 
(not shown), the region currently excluded by B s — > has my 2 < 180 GeV and 

m < 170 GeV, within the strips excluded by m^ and m x ± but allowed by b — > 37. We 
recall that the b — > 57 limit is very dependent on the assumption of universal scalar masses 
for the squarks, which does not play a role elsewhere in the analysis of constraints on CMSSM 
parameters, and is of course untested. Clearly B s — > fi + \i~ has the potential to complement 
or even, in the future, supplant the b — ► 37 constraint, though it also relies on squark-mass 
universality. 

The input value of the trilinear soft supersymmetry-breaking parameter Aq has a signifi- 
cant effect on the allowed CMSSM parameter space, and is also important for rrih as well as 
the b — *■ 57 and B s — > fi + \i~ decays. Therefore, we display in Fig. 3 the interplays of these 
constraints for tan/5 = 57, \i > 0, and (a) A = 2m 1 / 2 and (b) A = —2myz. The quali- 
tative conclusions are similar to the Aq = case discussed previously: the region currently 
disallowed by B s — > /i + fi~ decay largely overlaps with the regions previously disfavoured by 
m h and b — > 57, and decreases in extent as A is reduced from positive to negative values. 
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5 Treatment of Errors 



There are several ways to treat the auxiliary errors in the B s — > jj, + analysis. In the 
above, we have implicitly assumed one of the most conservative treatments, in the sense 
that it excludes the smallest region of the (mi/2,mo) plane for given fixed values of the 
uncertainties. However, other treatments are possible, and here we compare their results. 
For this comparison, we include all the uncertainties in f Bg , m t and m h discussed in the 
previous Section. 

In our previous treatment, see panel (d) of Fig. 2, we assumed that all the uncertainties 
have Gaussian error distributions, and defined the allowed region by discarding the upper 
2.5% tail of the likelihood distribution obtained by combining the experimental and auxiliary 
errors. This would be correct if the central value of the experimental measurement, after 
subtracting any backgrounds, was strictly zero. Alternatively, one might discard the upper 

5 % of the combined likelihood distribution. This would give the correct experimental upper 
limit if the central experimental value were far enough above zero that no significant part of 
the lower tail of the likelihood distribution extends below zero B s — > ji + ji~ decay rate, but is 
otherwise clearly more conservative than the previous treatment. Finally, experimentalists 
sometimes subtract one theoretical (systematic) error from the measured result and then 
plot the 95 % confidence-level contour given by the experimental error. 

The two alternative prescriptions yield similar upper limits on m , though the shapes of 
the allowed regions are different. They both reach up to m ~ 400 GeV, as compared the 
range m < 350 GeV found in the previous analysis, shown in panel (d) of Fig. 2. We prefer 
to use the more conservative prescription used in drawing the previous figures, also because 
it is demonstrably appropriate if the central experimental value is negligible, as is currently 
the case. 

6 Possible Future Developments 

It is expected that the Fermilab Tevatron collider experiments will continue to improve 
the present sensitivity to B s — > fi + decay. To assess the likely impact of this improved 
sensitivity, we exhibit in Fig. 4 the potential (m^, mo) planes for Aq = and tan j3 = 57, 40, 
obtained using the conservative error prescription and neglecting possible improvements in 
the determinations of fs s , m t and m&, in the pessimistic case that no signal is seen. 

In panel (a) of Fig. 4 for tan (3 = 57, we show that the B s — > fi + pT constraint would, 
under these pessimistic assumptions, begin to exclude a region in the neighbourhood of 
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Figure 4: The potential disallowed regions in the (mi/2,mo) plane for A = 0, \l > and 
(a) tan/? = 57, (b) tan/5 = 40, obtained assuming a Fermilab Tevatron upper limit on 
B s — > fi + that is improved to a 95% CL upper limit o/5x 10 -8 , with no parallel reductions 
in the uncertainties in Jb s , m t andrrib. Panels (c) and (d) show the corresponding disallowed 
domains assuming a conjectural LHC measurement BR(B S — > = (3.9 ± 1.3) x 10~ 9 . 

The various contours and shadings in the (mi/2,m ) plane are as explained in the text de- 
scribing Fig. 2. 
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(mi/2, m ) = (400, 400) GeV that is favoured by WMAP and allowed by all the other present 
constraints, assuming a future 95% CL upper limit of 5 x 10~ 8 . The region at small mi/ 2 
is allowed because of the high sensitivity to m t and m b seen in panel (b) of Fig. I. For 
tan/5 = 40, as seen in panel (b) of Fig. 4, the region disallowed by B s — > p + would 
still lie within the region already disallowed by b — > 57. The much reduced sensitivity of 
B s — > yU + decay to m t and m& at small mi/ 2 for tan (3 = 40 implies that there is no allowed 
'gap' at small m^. 

Panels (c) and (d) of Fig. 4 show the corresponding sensitivities at the LHC, assuming a 
conjectural 3-a measurement whose central value coincides with the central value predicted 
by the Standard Model, i.e., BR(B S -> fi+ fi~) = (3.9 ± 1.3) x 10~ 9 . In panel (c) for 
tan (3 = 57 we see that such a measurement would cover a very large fraction (but not all) of 
the CMSSM parameter space allowed by WMAP. On the other hand, the fraction covered 
in panel (d) for tan (3 = 40 is somewhat smaller. Of course, it is quite likely that the present 
errors in fs s , m t and m;, will be substantially reduced by the time of this ultimate LHC 
measurement. Reducing the fs 3 error, in particular, would reduce the scope available for a 
CMSSM contribution, and extend the B s — > fi + fi~ exclusion region to larger m^. 

7 Conclusions 

We have seen that the interpretation of the present and prospective experimental limits on 
B s — > fi + fx~ decay are very sensitive to auxiliary uncertainties, principally those in f Bs , m t 
and TOft. At the present time, these restrict significantly the regions of the CMSSM parameter 
space that can be excluded by the current upper limit on this decay, and their uncertainties 
may not be reduced significantly during the remaining operation of the Fermilab Tevatron 
collider, with the likely exception of m t . However, one might hope that the uncertainties 
in each of fs s , m t and m& could be reduced by the time the LHC achieves its ultimate 
sensitivity to B s — > /i + fj,~ decay. As an exercise, we have considered the possibility that their 
uncertainties might be reduced to Af Bs = 10 MeV, Am t = 1 GeV and Am b = 0.05 GeV. 
In this case, the LHC reach in the (mi/ 2 ,mo) plane for tan/5 = 57 would extend to m 1 / 2 ft 
1400 GeV along the WMAP strip. 

Beyond the framework of the CMSSM, it would be interesting to study the interpretation 
of the B s — > fi + n~ constraint also in the frameworks of more general models, such as those 
with non- universal Higgs masses (NUHM) 2 [42]. A complete study of the situation within 
the NUHM would take us beyond the scope of this paper, so we restrict ourselves to a few 
2 We expect the situation in the general low-energy effective supersymmetric theory (LEEST) to be similar. 
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remarks. The models likely to be disfavoured by B s — > fi + pT are those with a low value 
of vriA and large tan (5. Such models also tend to predict large neutralino-nucleus scattering 
cross sections [7]. We have examined a sample of NUHM scenarios that are apparently 
excluded by the recent CDMS II upper limit on the direct scattering of supersymmetric dark 
matter, and found that about half of them are excluded by B s — > fi + The next step 
would be to examine models apparently allowed by CDMS II, to see how many of them are 
also excluded by B s — > /i + However, for consistency, one should also study the effects on 
the dark matter scattering cross section of auxiliary uncertainties such as those in m t and 
m&, which has not yet been done in the manner described here for B s — > fi + or . We plan to 
present elsewhere such a unified treatment of the uncertainties. 
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